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In this paper we present the results of an optical and near infrared identification of 514 radio sources from the FIRST 
survey (Faint Images of the Radio Sky Survey at Twenty centimetres) with a flux-density limit of 1 mjy in the NOAO 
Deep- Wide Field Survey (NDWFS) Bootes field. Using optical (Bw, R, I) and K band data with approximate limits of 
Sw~ 25.5 mag, i?~25.8 mag, /~ 25.5 mag and A"~19.4 mag, optical counterparts have been identified for 378 of 
514 FIRST radio sources. This corresponds to an identification rate of 34% in four bands (BwRIK), 60% in optical 
bands (BwRl) and 74% in / band. Photometric redshifts for these sources have been computed using the hyperz 
code. The inclusion of quasar template spectra in hyperz is investigated. We note that the photometric data are, in 
many cases, best matched to templates with very short star-formation timescales and the inferred ages of identified 
galaxies depend strongly on the assumptions about the star-formation timescale. The redshifts obtained are fairly consis- 
tent with those expected from the K-z relation for brighter radio sources but there is more scatter in the K-z diagram at z< 1 . 
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1 Introduction 

Radio sources detected in deep 1 .4GHz surveys are gener- 
ally drawn from two populations (Condon, Cotton & Brod- 
erick 2002). Brighter sources are mostly active galactic 
nuclei (AGN) in which the radio emission is synchrotron 
emission generated by energetic charged particles interact- 
ing with magnetic fields. The particles are produced in the 
vicinity of the supermassive black holes which reside at the 
centre of galaxies. As one approaches 1 mjy flux-densities, 
source populations include starforming galaxies in which 
radio emission is a combination of synchrotron emission 
from charged particles, generated by star formation activ- 
ity, and thermal emission from hot gas in the environment 
of starforming regions. 

Understanding the sources found in 1 .4GHz surveys is 
thus relevant to many areas of study including star formation 
histories and processes, accretion and black hole physics 
and astrophysical jet formation. Moreover, radio sources 
can be detected out to very high redshifts over large ar- 
eas of sky and thus probe larger volumes than those typ- 
ically observed in optical surveys. This makes them valu- 
able cosmological probes. The clustering of these sources 
has been measured (Cress et al. 1996; Magliocchetti et al. 
1998; Overzier et al. 2002; Blake, Ferreira & Borrill 2004), 
as has the weak gravitational lensing (Chang et al. 2002) and 
these measurements can be used to put constraints on cos- 
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mology (Cress & Kamionkowski 1998; Negrello, Maglioc- 
chetti & De Zotti 2006; Kamionkowski et al. 1998), partic- 
ularly when the redshift distribution is well known. In addi- 
tion, if optical identifications allow sources to be separated 
into redshifts bins, then new information on the evolution of 
structure can be obtained. 

The largest previous sample of optical identifications 
for FIRST sources (Faint Images of the Radio Sky at 
Twenty centimeters; Becker, White & Helfand 1995) are 
based on the Cambridge Automated Plate-Measuring Ma- 
chine (APM) survey (Maddox et al 1990). An identification 
program described by McMahon et al. (2002) for 382, 892 
sources in the north Galactic cap, using the APM scans of 
POSS-I plates, was used as astrometric standard to improve 
the absolute astrometry of the POSS plates: matching the 
radio and the optical catalogues resulted in 70, 000 opti- 
cal counterparts with E< 20 and 0< 21.5 (where E and O 
bands are similar to R and B bands respectively). Maglioc- 
chetti et al. (2002) studied optical counterparts to FIRST 
sources in the 2-degree Field Galaxy Redshift Survey (2dF- 
GRS, Colless et al 2001). Magliocchetti & Maddox (2002) 
present a detailed analysis of the properties and angular 
clustering of ~ 4000 FIRST sources with bright optical 
counterparts in the APM survey. Ivezic at al. (2002) discuss 
the properties of ~ 30, 000 sources observed by the FIRST 
survey and the Sloan Digital Sky Survey (SDSS, York et al. 
2000). The ugriz observations of the SDSS have limiting 
AB magnitudes of 22.0, 22.2, 22.2, 21.3 and 20.5 respec- 
tively. SDSS Spectra are obtained for objects with r< 17.77 
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as well as for luminous red galaxies described in Eisenstein 
et al (2001). 

While the studies discussed above have provided large 
amounts of information on the optically bright counterparts 
to 1 .4GHz sources, one requires deeper surveys to identify 
the bulk of the radio sources. Optical counterparts to sources 
in the Phoenix Deep survey (Hopkins et al. 2003) have been 
studied by Sullivan et al (2004). They investigate photomet- 
ric redshifts for sources in a region with varying sensitivity 
- the theoretical rms noise ranges from 10 fiJy to 0.36mJy. 
They have U BVRI imaging for about 1 deg 2 but only have 
K-band observations for a 15' x 15' region in which they 
could find only a handful of sources with S > lmJy. Since IR 
observations provide important constraints for photometric 
redshift estimation at higher redshifts, little information is 
provided on the redshift distribution of the brighter sources. 
Similarly, in the work by Ciliegi et al. (2005) which inves- 
tigates optical counterparts to 80 fiiy sources, deep K-band 
data is only available in 165 arcmin 2 which allows poten- 
tially reliable photometric redshift estimation for a small 
number of sources and they do not address this explicitly. 

A complete subsample of the Leiden Berkeley Deep 
Survey (LBDS, Windhorst, van Heerde & Katgert 1984) 
has been studied with deep optical and IR observations. The 
LBDS Hercules subsample consists of 72 sources observed 
at 1.4 GHz with flux-density Si. 4 > lmJy in a 1.2 deg 2 
region of Hercules. They imaged sources in G,R,I and K 
bands, identifying 69 of the sources in at least two bands. 
They also obtained spectra of 47 of the sources (Waddington 
et al. 2000). Photometric redshifts were estimated for the 
rest of sample using the spectral synthesis code of Jimenez 
et al. (1998). The data was used to examine the evolu- 
tion of the radio luminosity function and its high redshift 
cut-off (Waddington et al. 2001). This study has allowed 
the redshift distribution and nature of 1 .4GHz sources with 
S> lmJy to be studied fairly well but the area of sky is small 
and may not represent a fair sample. The Combined EIS- 
NVSS Survey of Radio Sources (CENSORS) survey (Best 
et al. 2003, Brookes et al. 2006) matches sources in the ESO 
Imaging Survey to 150 sources in the NRAO VLA Sky Sur- 
vey (Condon et al. 1998) and is complete to 7.8mJy. Com- 
plete spectroscopic data from this survey will provide fur- 
ther information on the nature and redshift distribution of 
faint radio sources. 

Ultimately, we would like to study the evolution of clus- 
tering using measurements in different redshift slices and 
since it is very difficult to obtain spectra for a large number 
of optically faint sources, it is useful to investigate the use 
of photometric redshifts more extensively. Templates in the 
photometric redshift code used by Waddington et al. (2001) 
do not explicitly include light from AGN activity and it is 
useful to consider this (Brodwin et al. 2006). 

An interesting application of the optical identifications 
of radio sources involves using the K-z relation to probe 
the nature of radio galaxy hosts. Lilly and Longhair (1984) 
showed that for 3CRR galaxies, (sources with flux-densities 



greater than lOJy at 178MHz), there is a tight correlation 
between the K-band magnitudes of host galaxies and the 
redshift of the galaxy. The simplest interpretation of these 
results is that radio galaxies formed at high-redshift and 
have been evolving passively since then (see also Jarvis et 
al. 2001). Best et al. (1998) argued that this interpretation 
could not be correct since there is evidence that radio galax- 
ies at high redshift are found in richer environments than 
those at low redshift. Studies of fainter radio sources pro- 
vide more information on the nature of radio-galaxy hosts, 
with results from the 6C catalogues indicating an offset in 
the K-z relation to fainter magnitudes. Some studies (eg. 
Eales et al. 1997) claimed the offset only applies to z>0.6 
galaxies but Willott et al. (2003) attribute the apparent lack 
of offset at low-z to the small number of sources in the 
sample. Their results, as well as those of McLure et al. 
(2004), are consistent with the idea that the K-z relation 
for galaxies with brighter radio flux-densities (S>500 Jy at 
151MHz) can be understood using the correlation between 
stellar mass and black-hole mass (Magorrian et al. 1998) 
combined with a correlation between black hole mass and 
radio luminosity. Best et al (1998) and de Breuck et al. 
(2001) also use this idea in explaining their results. While 
the reasoning involving black hole mass correlations works 
well for brighter sources, in McLure et al.(2004), it appears 
to break down for the fainter radio sources in the TexOx- 
1000 survey (Hill & Rawlings 2003) and further studies of 
the K-z relation for fainter radio sources would help in un- 
derstanding this. 

Observations of extragalactic radio sources also provide 
important tests for galaxy evolution theory, particularly in 
clusters of galaxies where AGN are believed to heat the in- 
tracluster medium (eg Croton et al. 2006). 

In this paper, we discuss the properties and K-z rela- 
tion of optical counterparts of FIRST radio sources found 
in the Bootes field of the NOAO Deep-Wide Field Survey 
(NDWFS). Many studies are being carried out in this re- 
gion. Using FIRST matches to NDWFS sources, Wrobel et 
al. (2005) identified 55 candidates suitable for the study of 
parsec-scale properties of radio sources and McGreer et al. 
(2006) searched for high-redshift quasars, finding a z=6.1 
quasar. De Vries et al. (2002) and Wilman et al. (2003) stud- 
ied the properties of radio sources in the Bootes field de- 
tected in a Westerbork 1 .4 GHz survey and pointed out that 
optical and IR identifications were important for further in- 
terpretation of their results. New photometric data at Spitzer 
wavelengths (Brodwin et al. 2006) and X-ray wavelengths 
(Brand et al. 2006), as well as a spectroscopic campaign 
called the AGN and Galaxy Evolution Survey (Brown et al., 
2006) in the Bootes region, will provide further information 
on the radio sources in the field. The NDWFS data covers 
- 9 deg 2 (216.1° < RA < 219°, 32° < Dec < 36° ) but 
we only deal with the ~ 5deg 2 between 34° < Dec < 36° 
which has K-band data readily available. 

In section 2 and 3 we describe briefly the FIRST survey 
and the NDWFS. In section 4 we present the results from 
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the optical and infrared identification and, in section 5, we 
discuss the magnitude distributions of the optical counter- 
parts. Section 6 deals with the photometric redshift estima- 
tion, the resulting redshift distributions and the K-z relation. 
In section 7, we summarise our results. 

2 The FIRST Survey 

The FIRST Survey (Faint Images of the Radio Sky at 
Twenty centimetres) covers a total of 9033 deg 2 of the 
sky: 8422 deg 2 in the North Galactic cap and 611 deg 2 in 
the South Galactic cap. This yields a catalogue of 811117 
sources with a positional accuracy better than 1" at 90 % 
confidence and flux-density limit of 1 .0 mJy at 1 .4 GHz . We 
have used the 3 April 2003 version of the FIRST catalogue 
which is derived from the 1993 through 2002 observations. 
This catalogue is available from the FIRST website 1 . 




8 10 12 14 16 18 20 
Distance(arcsec) 

Fig. 1 Distribution of the radio-optical positional offsets 
between each FIRST radio source and the closest NDWFS 
object. The vertical dotted line indicates the distance cut-off 
(2 ) used to define the real radio-optical identifications. 



3 The NOAO Deep- Wide Field Survey 
Catalogue 

The NDWFS is a deep optical and near infrared (Bw RIJ H 
K) imaging survey that covers two 9.3 square degree fields 
designed to study the formation and the evolution of large 
scale structure (Jannuzi & Dey 1999). The first semester 
field (Bootes field ) is located near the north Galactic pole, 
centred on RA = U h 32 m 05.571 s , dec = +34° 16' 47.5" 
(J2000) ) and the second semester (Cetus field) is roughly 
30" from the South Galactic pole centred on RA of 2 h 10 m 
and dec = -4°30'. 

This paper utilises the NDWFS third data release (DR3) 
obtained for the 27 subfields that compromise the Bootes 
fields of the NOAO Deep Wide Field Survey. We have used 
the merged BwRIK catalogue for the Bootes field which is 
available in ASCII tables from the NOAO website 2 . The 
NDWFS object catalogues were generated using the SEx- 
tractor 2.3.2 (Brown et al. 2003). This catalogue has been 
split by declination range into four strips (32° < Dec < 
33°, 33° < Dec < 34°, 34° < Dec < 35°, 35° < Dec < 
36°) and band (BwRIK) into 16 files (each strip observed 
in four bands: Bw, R, I and K). The merged catalogue con- 
tains 2 SExtractor flags which provide information about 
the catalogue data quality: 1) FLAG_DUPLICATE =1 indi- 
cates that object is repeated elsewhere in the catalogue; 2) 
FLAG_SPLITMATCH= 1 is used for objects with uncertain 
matches (dubious matches or photometry). 

4 Optical and Infrared Identification of 
FIRST radio sources 

There are several techniques which could be used to cross- 
correlate radio and optical catalogues. One statistically ro- 

1 http://sundog.stsci.edu/first/catalogs/ 

2 http://www.noao.edu/noao/noaodeep/ 
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Fig. 2 Right ascension and declination distribution of the 
FIRST radio sources identified in all four bands in the ND- 
WFS Bootes field. The inset shows the coverage of the K- 
band data. 



bust method is the likelihood ratio technique of Sutherland 
& Saunders (1992). This method has been used often in or- 
der to identify radio sources (e.g. Gonzalez-Solares et al 
2004, Ciliegi et al 2005, Sullivan et al 2004, Afonso et al 
2005). However if the positional accuracies for both radio 
and optical catalogues are very high, positional coincidence 
alone can be adequate. Sullivan et al. (2004) demonstrated 
that the two methods produced very similar catalogues of 
matched objects in their study. We therefore positionally 
matched the FIRST radio sources with the two last strips 
of the NDWFS catalogue which cover ^5 deg 2 . We adopt 
a search radius of 2 arcsec as a good compromise when in- 
specting Fig 1, where the distance distribution of the radio- 
optical separation is presented. We use the K-band to esti- 
mate the contamination rate since they are the sparsest im- 
ages. The number of sources observed in the ^5 deg 2 of 
interest, in the K band is 9.2 x 10 4 (see table 1), giving a 
surface density, p, of 1.4 x 10~ 3 sources/arcsec 2 . The num- 
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1 st Strip 


2 nd Strip 


3 rd Strip 


4 th Strip 


Bw 


397531 


650832 


790446 


670230 


R 


339584 


606205 


605044 


565157 


I 


320125 


568260 


652376 


526129 


K 




3793 


51122 


40810 



Table 1 The number of sources in the NDWFS catalogue 
observed in Bw, R, I and K. Each column gives the number 
of sources in a 1° strip, the first strip being 32° < S < 33°, 
the second being 33° < 5 < 34° etc. Only the third and 
fourth strips are considered in this work. 



3 rd Strip 4 th Strip Both Strip 

34 < dec < 35 35 < dec < 36 34 < dec < 36 





Sources detected in one band 




Bw 


63% 


79% 


69% 


R 


64% 


85% 


72% 


I 


67% 


86% 


74% 


K 


37% 


46% 


40% 




Sources detected in two bands 




BwR 


58% 


74% 


63% 


Bwl 


57% 


72% 


63% 


RI 


60% 


78% 


67% 


IK 


35% 


44% 


38% 




Sources detected in three bands 




BwRI 


54% 


69% 


60% 


BwIK 


32% 


39% 


35% 


RIK 


34% 


43% 


38% 




Sources detected in four bands 




BwRIK 


32% 


38% 


34% 



Table 2 Optical and near infrared identification results for 
the FIRST radio sources. 



ber of FIRST sources in the region, Np, is 514. Thus, with 
a search radius of r s , the number of random coincidences is 
given by 

N r = N f x 2tt x r 2 x p w 18. (1) 

This gives a contamination rate of 3.5%. As a compari- 
son, Magliocchetti & Maddox (2002) used a 2 arcsec cutoff 
for the APM-FIRST matches corresponding to 5% random 
matches and Prandoni et al. (2001) used a 3 arcsec cutoff for 
the ATESP-EIS (Australia Telescope ESO Slice Project and 
ESO Imaging Survey (EIS)) with 2% contamination rate. 

We removed sources with uncertain matches in the 
merged catalogue (FLAG_SPLITMATCH= 1 ) . We also used 
the SExtractor star-galaxy classifier to investigate subsam- 
ples of the matches: we considered a source to be extended 
if the stellarity parameter was < 0.5. 

Figure 2 shows the position of FIRST sources identified 
in all four bands of the NDWFS Bootes field. 

Table 1 shows the number of NDWFS sources in each 
band in each strip. Table 2 shows the results of the opti- 
cal and IR identifications for the third and fourth strips in 
one band, two bands, three bands and four bands. The ar- 
eas covered by the third and fourth strips are approximately 



equal yet there are 198 FIRST sources in the third strip and 
316 FIRST sources in the fourth strip. The third strip also 
contains ~15% more sources in all bands. Interestingly, the 
identification rate is higher in the fourth strip where the 
source density is lower. This may be explained by the pres- 
ence of large-scale structure in the region of the third strip, 
where many of the galaxies detected in the radio are too 
faint to detect optically - either because they are obscured 
by dust or are too distant - but systematics in the surveys 
should also be considered. 

5 Magnitude distribution 

In Figure 3, we compare magnitude distributions of sources 
identified in all four bands. Moving from top to bottom, 
the histograms shown are distributions for (i) all NDWFS 
sources in the region in the given band, (ii) the FIRST 
sources identified in the given band and (iii) the FIRST 
sources identified in all four bands. Comparing the first and 
second histograms, one notes that the number of FIRST 
galaxies per magnitude interval does not continue to in- 
crease down to the magnitude limit, but rather turns over 
at Bw ~ 25, r~ 24, I ~ 23 and K~ 18. The turnover 
is <~ 0.5 magnitudes brighter than the magnitude limit for 
Bw, about a ~ 1.5 magnitudes brighter for R and ~ 2.5mag 
brighter for I and K, highlighting the redness of the pop- 
ulation of sources associated with mJy radio sources. The 
K-band is clearly shallower than the other bands. 

In Figure 4, we investigate sources which are resolved 
in the NDWFS (have a SExtractor stellarity parameter of 
less than 0.5). Moving from top to bottom, the histograms 
show the distributions for (i) all NDWFS sources in the re- 
gion in the given band, (ii) for all the NDWFS sources in 
the region with galaxy identifications in the given band, (iii) 
the FIRST sources identified in the given band and (iv) the 
FIRST sources identified in all four bands. Comparing the 
top two histograms, one sees that fainter sources are less 
likely to be characterised as galaxies, presumably because 
it is difficult to resolve more distant and intrinsically fainter 
sources, otherwise the trends are similar to those seem in 
Figure 3. 

It is interesting to note that in Waddington et al. (2001) 
54 out of sources 72 (75%) have K-band matches with 
K<19.4 while only 41% of sourcs are idenitfied here. This 
supports the claim that studying areas larger than <~1 deg 2 
is required to get a fair sample of mJy sources and is per- 
haps to be expected when one considers that radio sources 
are clustered. 

6 Photometric redshift distribution 

We estimated photometric redshifts for radio sources iden- 
tified in all the four bands: Bw, R, I and K, using the hyperz 
code developed by Bolzonella, Miralles & Pello (2000). The 
code takes a series of galaxy spectrum templates, varies the 
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20 22 24 26 28 
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Fig. 4 Magnitude distributions in the Bw, R, I and K bands for, from top to bottom, all sources in the NDWFS, sources 
in the NDWFS which are extended (have stellarity index < 0.5), the extended optical counterparts of the FIRST radio 
sources identified in the one band, the extended optical counterparts of the FIRST sources which have identifications in all 
four bands. 
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2 3 4 
Redshift z 




2 3 4 
Redshift z 




2 3 4 
Redshift z 




2 3 4 
Redshift z 



Fig. 5 Photometric redshift histograms for sources with counterparts in the four bands, as calculated by hyperz. Top 
left: Quasar templates and burst templates included, M>-27. Top right: Quasar templates included, no burst templates, 
M>-27. Bottom left: No quasar templates, burst templates included, M>-24. Bottom Right: No quasar templates, no burst 
templates, M>-24. 



redshift and amount of reddening and calculates photomet- 
ric measurements that would be made through selected fil- 
ters. These are matched to an input catalogue of photomet- 
ric measurements using a standard \ 2 minimisation proce- 
dure. We used GISSEL 98 galaxy templates which are the 
1998 update of the spectral synthesis models described by 
Bruzual & Chariot (1993). All the models have solar metal- 
licity and Miller & Scalo (1979) initial mass function. Each 
model corresponds to a different present-day spectrophoto- 
metric galaxy type. Galaxies are represented by an exponen- 
tially decaying star formation rate with a time scale of 1 Gyr 
for elliptical galaxies, increasing to 30 Gyr for Sd galaxies. 
Irregular galaxies are represented by a constant star forma- 
tion rate. The code also offers 'delta-burst' templates cor- 
responding to a single burst of star formation of a specific 
age. 

We adopted the reddening law of Calzetti et al. (2000). 
While there are large number of filters predefined in hyperz, 
the Bw filter is not included, we added the profile transmis- 
sion of the filter from the NDWFS website. We also inves- 
tigated setting a maximum absolute magnitude for sources. 
The default limit is -29 while brightest cluster galaxies have 
a maximum V-magnitude of about -23.5. We limited the ab- 
solute magnitudes to -24, except when including QSO spec- 
tra when we allowed the absolute magnitude to be as large 
as -27. 

We investigated the addition of four QSO templates 
obtained from the LE PHARE website 3 . These template 

3 http://www.lam.oamp.fr/arnouts/LE_PHARE.html 



where selected for only 9/177 sources (5%) and 5 of these 
were the highly reddened quasar template. 131 sources 
(74%) where identified as 'bursts', 14 sources (8%) where 
identified with ellipticals and 8 (5%) with SO's leaving only 
15 sources (8%) identified with late-type galaxies. When the 
'burst' option was removed, the reddest QSO template was 
selected as the best fit for 57/177 sources (32%) and the el- 
liptical template was selected for 76 (43%) of sources. Prob- 
abilities for QSO matches and elliptical matches were sig- 
nificantly lower than those for burst options in most cases. 

When we limited the absolute magnitudes to -24 and 
removed the QSO templates. 138 of the sources were iden- 
tified with 'bursts', 15 with ellipticals, 9 with SO's and 15 
with later-type galaxies. When burst spectra and QSO spec- 
tra were not available, 138 sources (78%) were identified 
with ellipticals. 

As in Bolzonella et al. (2000) we interpret the 'burst' 
spectra as elliptical galaxies with star formation on 
timescales shorter than 1 Gyr. When sources are identified 
with bursts, their ages range from 0.003 to 9.5 Gyr with 
a mean of 1.3 Gyr. When the burst option is removed and 
sources are identified with ellipticals, ages range from 0.01 
to 12.5 Gyr with a mean age of 5.6 Gyr It is interesting to 
note that the age-to-formation implied for these 'ellipticals' 
is much longer than the ages implied for 'burst' spectra, sug- 
gesting that strong conclusions based on photometric age 
dating of ellipticals should be treated with caution. 

We also considered the subset of identifications which 
were resolved (SExtractor stellarity parameter < 0.5). 
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Fig. 6 The K-band magnitude versus redshift for FIRST- 
NDWFS matches. The upper line gives our fit to the FIRST 
data, the more curved line is the fit to 7C sources, similar 
to that obtained by Brookes et al (2006) for the CENSORS 
survey, and the third line is the fit to the 3CRR, 6C and 7C 
combination Willott et al. (2003). 



When QSO spectra are offered, 4% are identified as such. 
When we limit the absolute magnitude to -24, 94/1 13 (83%) 
of sources are then identified with the 'delta burst' templates 
with a mean age of 1 .6 Gyr. Four sources are identified with 
ellipticals, 5 with SO's and 10 with later type galaxies. The 
ages and split between different types of galaxies are similar 
to that obtained for all sources. 



We tried to estimate photometric redshifts for sources 
that were only detected in three bands but found, by drop- 
ping one measurement for sources that were detected in all 
four bands, that estimates using three bands were not reli- 
able. 



Figure 5 shows the photometric redshift distribution for 
the FIRST radio sources identified in the NDWFS Bootes 
field when various templates are included. The magnitude 
cut-off is -24 when quasars are not included and -27 when 
they are included. If we allowed quasar templates to be 
matched with a cutoff of -27 but require other templates 
to be matched with a -24 cutoff, the distribution looks very 
similar to that shown in the bottom left panel (since there are 
very few sources identified as quasars). One notes, in the top 
right panel, that matching to the red quasar template when 
bursts are not available, results in increase in sources in the 
l<z<1.5 range. The burst option also tends to increase the 
number of lower redshift objects. 



7 The K-z relation 

In Figure 6 we plot the K-band magnitude of FIRST sources 
against photometric redshifts obtained for sources identified 
in all four bands of NDWFS (including unresolved sources). 
The redshifts correspond to those given in the bottom left 
panel of Figure 5 where a magnitude cutoff of -24 is used. 

The K-z relation for the combination of radio galax- 
ies in 3CRR, 6C and 7C catalogues, given in Willott et 
al. (2003), is plotted in the figure (K=17.37 + 4.531og 10 z- 
0.31 (log 10 z) 2 ). The fit for 7C alone, which Brookes et al 
(2006) claim is similar to the fit they obtain for the CEN- 
SORS survey, is also plotted. A fit to our data, given by 
K=(17.9±0.1)+(4.3±0.3)logi z, is the line at faintest mag- 
nitudes. We have excluded the clump of sources at z=0.05 
in the fit as these sources appear to be spurious identifica- 
tions, possibly resulting from a lack of resolution in redshift 
during template matching. They are mostly identified with 
starforming galaxies and would not be expected to fit the 
K-z relation for radio galaxies. In the fit, we do not use the 
uncertainties on the redshifts generated by hyperz as these 
are found to be unreliable when tested on samples where 
spectroscopic redshifts are available. Since FIRST sources 
are much fainter than 7C sources, which have flux-densities 
Si5iMHz > 500 mJy, and are selected at a much higher 
frequency, one might not expect the established K-z rela- 
tion to hold. In addition, the uncertainties in the photomet- 
ric redshift estimates are large and different aperture sizes 
for photometry could affect the results, so it is interesting 
that our relation is so similar to the other two. The fit for 
the FIRST sources is shifted slightly to fainter magnitudes 
which could support findings that brighter radio sources are 
associated with galaxies that are brighter in K, even for faint 
radio sources, but a more complete sample is required to in- 
vestigate this further. We considered FIRST sources with 
radio flux-densities greater than lOmJy and did not find a 
better fit to the Willott et al. (2003) relation. 

In Figure 7, we plot the magnitude in each band against 
photometric redshift and find signs of a similar correla- 
tion in other bands, except with a larger scatter, particu- 
larly in I. The la variations are given by gb = 2.65, gr = 
2.37, 07 = 9.77, <j K = 1-96. It is unclear why the I-band 
should have particularly large scatter. We note that the scat- 
ter in the K-band is significantly smaller at z>l than at 
lower redshift, consistent with the idea that fainter radio 
surveys include fainter galaxies that are different from the 
typical galaxies associated with bright radio sources. Once 
again, a more complete sample and more reliable estimates 
of redshifts, which include more spectroscopic measure- 
ments, are required to investigate this further. 

8 Conclusion 

In this paper we have presented the optical and near in- 
frared identifications of 514 radio sources in FIRST survey 
obtained by matching objects in the NDWFS survey, over 
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Fig. 7 Bw magnitude versus redshift (top left), R magnitude versus redshift (top right), I magnitude versus redshift 
(bottom right) and K magnitude versus redshift (bottom left) for the identified radio sources. 



the region 216.1° < RA < 220°, 34° < DEC < 36°. We 
identified 177 sources (35%) in all 4 bands of NDWFS. In 
the I-band, 74% of sources were identified. We found a sur- 
prisingly large difference between the number of sources in 
the upper declination strip and that in the lower strip, possi- 
bly explained by the presence of large-scale structure in the 
34° < DEC < 35° strip. Spectroscopy would be required to 
explore this further. Magnitude distributions demonstrated 
the red colours of the identified galaxies. The subset of 
sources which are extended (stellarity <0.5) in all bands 
(113/177) showed similar trends to the set including point 
sources but missing, perhaps, contributions from high-z ob- 
jects which are difficult to resolve. 

We computed photometric redshifts for those sources 
with identifications in four bands using the code hy- 
perz. The 'burst' templates, corresponding to populations 
which formed over very short timescales, provided the best 
matches for the majority of sources. The ages implied 
for galaxies which match burst templates are significantly 
shorter than those implied if they are matched to elliptical 
templates. The four quasar templates we tried did not ap- 
pear to fit the colours very well, although the dust-reddened 
QSO template was selected as the best match in many cases 
when the burst templates were not available. It appears that 
alternative approaches to photometric redshift estimation, 
such as those described in Brodwin et al. (2006) need to be 
explored. 

The K-z relation we obtained using photometric red- 
shifts is similar to that obtained for brighter sources inves- 
tigated by Willott et al. (2003) and Brookes et al. (2006) 
although there is much more scatter below z=l. Deeper K- 



band data and more use of spectroscopic redshifts is re- 
quired before strong conclusions can be drawn. We plan to 
carry out a similar analysis for sources in the southern ND- 
WFS field, which we can follow up spectroscopically using 
the Southern African Large Telescope. 
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